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Abstract —In actual implementations of magnetic control laws 
for spacecraft attitude stabilization, the time in which Earth 
magnetic field is measured must be separated from the time in 
which magnetic dipole moment is generated. The latter separation 
translates into the constraint of being able to genere only 
piecewise-constant magnetic dipole moment. In this work we 
present attitude stabilization laws using only magnetic actuators 
that take into account of the latter aspect. Both a state feedback 
and an output feedback are presented, and it is shown that the 
proposed design allows for a systematic selection of the sampling 
period. 

Keywords—attitude control of spacecraft; magnetic actuators; 
control algorithms implementation; averaging; Lyapunov stability. 

I. Introduction 

Spacecraft attitude control can be obtained by adopting 
several actuation mechanisms. Among them electromagnetic 
actuators are widely used for generation of attitude control 
torques on satellites flying low Earth orbits. They consist of 
planar current-driven coils rigidly placed on the spacecraft 
typically along three orthogonal axes, and they operate on 
the basis of the interaction between the magnetic moment 
generated by those coils and the Earth’s magnetic field; in fact, 
the interaction with the Earth’s field generates a torque that 
attempts to align the total magnetic moment in the direction 
of the field. Magnetic actuators, also known as magnetorquers, 
have the important limitation that control torque is constrained 
to belong to the plane orthogonal to the Earth’s magnetic 
field. As a result, sometimes a reaction wheel accompanies 
magnetorquers to provide full three-axis control (see HI Sec¬ 
tion 7.4]); moreover, magnetorquers are often used for angular 
momentum dumping when reaction or momentum-bias wheels 
are employed for accurate attitude control (see IT] Section 
7.5]). Lately, attitude stabilization using only magnetorquers 
has been considered as a feasible option especially for low- 
cost micro-satellites and for satellites with a failure in the 
main attitude control system. In such scenario many control 
laws have been designed, and a survey of various approaches 
adopted can be found in ||2|; in particular, Lyapunov based 
design has been adopted in ||3l, ID, Q, ||6l, Q. In those works 
feedback control laws that require measures of both attitude 
and attitude-rate (i.e. state feedback control laws) are proposed; 
moreover, in 0 and m feedback control algorithms which 
need measures of attitude only (i.e. output feedback control 
algorithms) are presented, too. 


ment of the geomagnetic field and a continuos-time generation 
of coils’s currents. However, in practical implementations the 
time in which Earth’s magnetic field is measured (measurement 
time) has to be separated from the time in which coils’ 
currents are generated (actuation time). The latter separa¬ 
tion is necessary because currents flowing in the spacecraft’s 
magnetic coils generate a local magnetic field that impairs 
the measurement of the geomagnetic field; as a result, when 
the Earth magnetic field is being measured no currents must 
flow in the coils, and consequently no magnetic actuation is 
possible. As a result, in practical applications a periodic switch 
between measurement time and actuation time is implemented. 
Usually measurement time is set much shorter than actuation 
time, so that the spacecraft is left unactuated for a very short 
fraction of the period; thus, the measurement process can be 
modeled as instantaneous, and during each period the magnetic 
dipole moment is kept constant since only the measurement at 
the beginning of the period is available. Thus the separation 
illustrated before translates into the constraint of being able to 
generate only constant magnetic dipole moments during each 
period. Then, magnetic control laws compatible with the latter 
constraint could be obtained by simply discretizing continuos- 
time control algorithms (see e.g. 0); however, here a different 
approach is proposed; it consists in taking into account of the 
latter constraint directly during the control design phase. It 
will be shown that the proposed approach leads to a systematic 
selection of the sampling interval; the latter aspect represents 
an advantage with respect to discretization of continuous-time 
control laws where the sampling interval is usually determined 
by trial an error. The proposed method is presented for both a 
state and an output feedback . 


A. Notations 


Eor X G K", ||a;|| denotes the Eucledian norm of x; for 
matrix A, ||A|| denotes the 2-norm of A. Symbol / represents 
the identity matrix. Eor a G represents the skew 

symmetric matrix 
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so that for b G multiplication a^b is equal to the cross 
product a X b. 


II. Modeling and control algorithms 


All control laws presented in the cited works are continuos- 
time; thus, in principle they require a continuous-time measure¬ 


In order to describe the attitude dynamics of an Earth¬ 
orbiting rigid spacecraft, and in order to represent the geo- 




magnetic field, it is useful to introduce the following reference 
frames. 

1) Earth-centered inertial frame Ti. A commonly used 
inertial frame for Earth orbits is the Geocentric Equa¬ 
torial Erame, whose origin is in the Earth’s center, its 
Xi axis is the vernal equinox direction, its Zi axis 
coincides with the Earth’s axis of rotation and points 
northward, and its yi axis completes an orthogonal 
right-handed frame (see |[T] Section 2.6.1] ). 

2) Spacecraft body frame Fb- The origin of this right- 
handed orthogonal frame attached to the spacecraft, 
coincides with the satellite’s center of mass; its axes 
are chosen so that the inertial pointing objective is 
having Fb aligned with Fi. 


Since the inertial pointing objective consists in aligning 
Fb to Fi, the focus will be on the relative kinematics and 
dynamics of the satellite with respect to the inertial frame. 
Let q = [qi 92 Qs = [qI QaV with ||g|| = 1 be the unit 
quaternion representing rotation of Fb with respect to Fi ; then, 
the corresponding attitude matrix is given by 

C{q) = [ql - qlqv)l + ‘ 2 qvql - 2949 ^ ( 2 ) 

(see 0 Section 5.4]). 


Let 


W{q) 


1 

2 


qd + qf 


(3) 


Then the relative attitude kinematics is given by 

q = W{q)uj (4) 

where a; G is the angular rate of Fb with respect to Fi 
resolved in Fb (see m Section 5.5.3]). 


The attitude dynamics in body frame can be expressed by 

Juj =—uj^ Juj + T (5) 

where J G is the spacecraft inertia matrix, and Q G 

is the control torque expressed in Fb (see ||9] Section 6.4]). 

The spacecraft is equipped with three magnetic coils 
aligned with the Ft axes which generate the magnetic control 
torque 

Q — trtcoiis X E? — B rriQQiis (6) 

where rricoiis G is the vector of magnetic dipole moments 
for the three coils, and is the geomagnetic held at spacecraft 
expressed in body frame Fb- 

Let be the geomagnetic held at spacecraft expressed 
in inertial frame Fi. Note that varies with time at least 
because of the spacecraft’s motion along the orbit. Then 

B\q,t) = C{q)B\t) (7) 

which shows explicitly the dependence of B^ on both q and 
t. 

Grouping together equations © 0 (§ the following 
nonlinear time-varying system is obtained 

q = W{q)uj 

JUJ = -CO^JUJ - B’’{q,ty TTlcoils ^ 

in which rricoiis is the control input. 


In order to design control algorithms, it is important to 
characterize the time-dependence of B^{q, t) which is the same 
as characterizing the time-dependence of B'‘{t). Assume that 
the orbit is circular of radius i?; then, adopting the so called 
dipole model of the geomagnetic held (see ifTOl Appendix H]) 
we obtain 

B\t) = ^Mm^yr\t))r\t)-m^] (9) 

In equation ©, pm = 7.746 10^^ Wb m is the total dipole 
strength (see Ini '). F(t) is the spacecraft’s position vector 
resolved in Fi, and F{t) is the vector of the direction cosines 
of r*(f); hnally is the vector of direction cosines of the 
Earth’s magnetic dipole moment expressed in Fi. Here we set 
TO* = [0 0 — 1]^ which corresponds to setting the dipole’s 
colelevation equal to 180°. Even if a more accurate value for 
that angle would be 170° (see HD), here we approximate it 
to 180° since this will substantially simplify future symbolic 
computations. 

Equation (j9]l shows that in order to characterize the time 
dependence ofi?* (t) it is necessary to determine an expression 
for r*(f) which is the spacecraft’s position vector resolved in 
Fi. Dehne a coordinate system Xp, Pp in the orbital’s plane 
whose origin is at Earth’s center; then, the position of satellite’s 
center of mass is clearly given by 

xP{t) = Rcos{ntfo) 

y^{t) = Rsm{nt + fo) ^ 

where n is the orbital rate, and fo an initial phase. Then, 
coordinates of the satellite in inertial frame Fi can be easily 
obtained from ( [T0| ) using an appropriate rotation matrix which 
depends on the orbit’s inclination incl and on the right 
ascension of the ascending node H (see fl] Section 2.6.2]). 
Plugging into (j^ the expression of the latter coordinates, an 
explicit expression for B'^{t) can be obtained; it can be easily 
checked that B'^{t) is periodic with period 271 jn. Consequently 
system © is a periodic nonlinear system with period 27: jn. 

As stated before, the control objective is driving the space¬ 
craft so that Fb is aligned with Fi. Erom © it follows that 
C{q) = I for q = [q^ qf\^ = ±q where q 0 0 1]^. Thus, 
the objective is designing control strategies for rricoiis so that 

—>■ 0 and w —>• 0. 

In Q both a static state feedback and a dynamic output 
feedback are presented; both feedbacks were obtained as 
modifications of those in a, 0 . 

The static state feedback control proposed in [|7] is given 
by 

rricoiisit) = {B'^{q{f),t)'^Y{e^kiqy{f) -h efc 2 w(f)) (11) 

and it is shown that picking ki > Q, k^ > 0, and choosing 
e > 0 and small enough, local exponential stability of {q, ui) = 
(q, 0) is achieved. 

The dynamic output feedback proposed in Q is given by 

S(t) = a(q(t) — eAS(t)) 

rricoiisit) = iB'’{qit),t)^)'^e‘^ {kiqy{t) ( 12 ) 

-\-k 2 aXW{q{t))'^{q{t) - e\S{t))) 

with i5 G K"*. It is shown that selecting fci > 0, /c 2 > 0, a > 0, 
A > 0, and choosing e > 0 small enough, local exponentially 
stability of equilibrium {q,u},d) = {q,0, -^q) is achieved. 





Both control laws (fTTli and (12 1 seem implementable in 
practice since ijcan be measured using magnetome¬ 

ters, qv{t) can be measured using attitude sensors, and uj{t) can 
be measured by attitude rate sensors. However, as explained in 
the introduction, the time in which B’^ is measured should be 
separated from the time in which control action is applied; in 
fact, magnetic control torque is obtained by generating currents 
flowing through magnetic coils, and those currents create a 
local magnetic field which impairs the measurement of B^. 
As a result, in order to take into account of the previous 
constraint, we should consider only feedback laws in which 
B’’ is instantaneously measured at the beginning of some fixed 
length interval, rricoiis is held constant over that interval, and 
the latter operations are repeated periodically. In the sequel T 
will denote the sampling interval, that is the length of each of 
those intervals. 


A static state feedback law that fulfills the previous con¬ 
straint can be designed by simply discretizing fn) through 
a sample-and-hold operation, thus obtaining the following 
control algorithm 


mcoiis{t) = {B^iqikT),kT)^)'^{e'^kiqy{kT) + ek2Uj{kT)) 
fc = 0,1,2,... kT<t<{k + l)T (13) 


Similarly, a dynamic output feedback law compatible with 
the constraint of generating a piecewise-constant signal for 
rricoiis can be obtained by discretizing ( [T^ using the forward 
differencing approximation (see ifT^ l. thus obtaining the fol¬ 
lowing control law 

6{{k + l)T) = S{kT)+Taiq{kT)-eXS{kT)) 
mcoUsit) = {B\q{kT),kTYye^{kiqy{kT) 

+k 2 a\W{q{kT)f{q{kT) - eX6{kT))) 
k = 0,1,2,... kT<t<ik + l)T 

(14) 

with 5 S 

III. State-feedback Design 

In this section we will focus on state-feedback law 
and show how to choose the parameters ki, k 2 , e, and the 
sampling interval T so that equilibrium {q, uj) = {q, 0) is 
locally exponentially stable for closed-loop system 
In order to derive those design criteria, it suffices considering 
the restriction of closed-loop system to the open set 

X where 

S3+ = {g G I ii^ii = > 0} (15) 


where 


Wy{qy) = ^[{l-qlqyy^^I- 


(17) 


and 

Rviqv) = (1 - 2qlqy) I + 2g„gT _ 2 (l - q^qy)^^^ q^ (18) 

and where 0 has been used. It is simple to show that if 
(g^,w) = (0,0) is a locally exponentially stable equilibrium 
for ( [T^ , then (g, w) = (g, 0) is a locally exponentially stable 
equilibrium for ([^ ( [T3] l. 

Next, consider the linear approximation o f ([T6] l about the 
equilibrium {qy,uj) = (0,0) as defined in 1131 . The latter 
approximation is given by 


Qv {t) 
Juj{t) 


-B'‘{t)^mcoiisikT) 


kT < t < {k + 1)T 
(19) 


rricoasikT) = B\kT)^^ {e^k^qyikT) + ek2u;{kT)) 

fc = 0,1, 2,... (20) 


Note that since is bounded (see 0), then Theorem II. 1 in 
applies to the nonlinear time-varying hybrid system 
consequently, {qy,uj) = (0,0) is a locally exponentially stable 
equilibrium for (16 if and only if it is an exponentially stable 
equilibrium for (19 1 ([20]l. 


Next consider the continuous-time system (T^ and sample 
its state [gt,(t)^a;(f)^]^ at f = (A:-|-1)T. The following discrete¬ 
time system is thus obtained 


g^((fc-|-l)r) = qy{kT) +'^uj{kT) 

- J~^ Gi{k,T)mcoiis{kT) 

uj{{k + l)T) = w{kT) - J-^G2{k,T)mcoUs{kT) 

k = 0,1,2... (21) 


where 


j■(k+l)T , 

Gr(k,T)= ((fc + l)T-r)S®(r)Xdr 

JkT ^ 


and 


Ak+l)T 

G2{k,T)= / B^rYdr 

JkT 


By ifTrl Proposition 7] it follows that if the linear time-varying 
discrete-time system ( [20| pT| ) is exponentially stable then the 
linear time-varying hybrid system ([Tg @ is exponentially 
stable. 


Since on the latter set q^ = (1 — gjg„)^^^ then the considered 
restriction is given by the following nonlinear time-varying 
hybrid system 

qv{t) = Wy{qy{t))Uj{t) 

JCoit) = -[Ry{qy{t))B\t)Ymcoiis{kT) 

kT<t< {k + l)T 

rricoUsikT) = [Ry{qyikT))B\kT)]^^ (e^kiqyikT) 

+ ek2Uj{kT)) 

fc = 0,l,2,... (16) 


Based on the previous discussion, our objective has become 
studying stability of linear time-varying discrete-time system 
( [20| ( [2T| . For that purpose it is useful to perform the following 
change of variables zi{k) = qy{kT), Z 2 {k) = uj{kT)/e 
obtaining 


zi{k + 1) 
Z2{k + 1) 


zi(fc) + ef Z 2 (fc) - eV-iGi(fc,T)B®(fcr)x^ 
(fclZi(fc) -I- k2Z2{k)) 

Z2{k) - eJ-^G2{k, T)B\kTY^ 

(fclZi(fc) -I- k2Z2{k)) 

( 22 ) 









Note that since Gi{k,0) = 0 i = 1,2, then Gi{k,T) can be 
expressed as 

G,ik,T) = H,ik,T)T 1 = 1,2 (23) 

where 

/■i 

H,{k,T)= -^{k,sT)ds 1 = 1,2 

Let L{k,T) = H 2 {k,T)B^{kT)^^. Note that since B^{t) is 
periodic, then the following average 

^ k—s-\-N 

Lav{T)=\iiR — L{k,T) 

k—s+1 


then 


Vi{wi,W2) = -k2wlL° 


,W2 


By La Salle’s invariance principle ifT^ Theorem 4.4], it follows 
that (27 1 is exponentially stable, and thus is Hurwitz. Then, 


by continuity there exists T* > 0 such that As{T) is Hurwitz 
for all 0 < r < T*. Next, fix 0 < T < T*, and consider 
symmetric matrix Pg > 0 such that 


PsAg{T) + Ag{TyPg = -I 


Let V 2 {z) = z^PgZ be a candidate Lyapunov function for 
system ([25]l. Then, the following holds 


is well defined and indipendent of integer s. The expression 
of Lav has been computed symbolically using Mathematical*^, 
but here it is omitted to save space. It turns out that = 
limT-i-o Lav{T) is symmetric, and that is positive definite 
if and only if the orbit is not equatorial, that is if its inclination 
satisfies incl ^ 0. Thus, we make the following assumption. 


Assumption 1: The spacecraft’s orbit satisfies condition 

Lt, > 0. 


At this point, consider the average system of (22 1 as defined 
in uni which is given by 

zi{k + l) = zi(fc) + ef Z 2 (fc) 

Z2{k + 1) = Z2{k) - (TJ~^Lav{T){kiZi{k)+ k2Z2{k)) 

(24) 

It can be verified that all assumptions of ifTSl Theorem 2.2.2] 
are fulfilled by system (22 1 . Thus, the following proposition 
holds true. 


Proposition 1: If the average system ( [24) i is exponentially 
stable for all 0 < e < eo, then there exists 0 < ei < eo^ such 
that system (22i is exponentially stable for all 0 < e < ei. 

For the average system ( |24] l the following stability result 
holds. 


Theorem 1: Assume that fci > 0, ^2 > 0 and Assumption 
is satisfied. Then there exists T* > 0, and for every 0 < 
T < T*, there exists eg > 0, such that fixed 0 < T < T* 
system ( p4| ) is exponentially stable for all 0 < e < cq. 

Proof: Rewrite ( |24l i in the following compact form 

z{k + l) = z{k) + eTAg{T)z{k) (25) 


with 

MT) 


0 

-kiJ-^Lav{T) 


-k2J-^Lav{T) 


(26) 


It will be shown now that = limT->.o Ag (T) is a Hur¬ 
witz matrix. For that purpose it is useful to introduce the 
continuous-time system w = A^w which in expanded form 
reads as follows 


1 

Wl = -W2 

JW2 = -L*^a{klWi + k2W2) 


(27) 


Introduce the Lyapunov function for ( [27| 

Vi{wi,W2) = kiwjLlaWi + ^wIJW2 


AV 2 {z) = (z + tTAg{T)zfPg{z + eTAg{T)z) - z'^PgZ 
< eT{-l + eT\\AgiTfPgAgiT)\\)\\z\\^ 


As a result, setting 


1 

2T\\Ag{T)-^PgAg{T)\\ 


(28) 


we obtain that system ( |25] l is exponentially stable for all 0 < 
e < eo. 


Based on the previous theorem, fix 0 < T < T*, and 
consider the corresponding value of eg given by ( [28] l; from 
Proposition 2 it follows that there exists 0 < ei < eo such 
that system ( [22| is exponentially stable for all 0 < e < ei. In 
conclusion, from Theorem [T] Proposition [T] and the preceding 
discussion, the following main proposition is obtained. 

Proposition 2: (State feedback). Consider the magnetically 
actuated spacecraft described by and apply the piecewise- 
constant static state-feedback ([13 i with fci > 0 and /c 2 > 0. 
Then, under Assumption [T] there exists T* > 0, and for every 
0 < T < T*, there exists ei < eg, with eg given by ( |28] l, 
such that fixed 0 < T < T* , for all 0 < e < ci equilibrium 
{q, w) = {q, 0) is locally exponentially stable for (|^ ( [T3| l. 

IV. Output-feedback Design 

In this section we will focus on output-feedback law ( [l4| 
and show how to choose parameters fci, k 2 , a. A, e, and 
sampling interval T so that equilibrium {q, uj, 5) = {q, 0, 
is locally exponentially stable for closed-loop system (j^ ( |14| l. 

Following the same approach as in state-feedback design, 
we determine the linear time-varying discrete-time system 
given by the interconnection of © with the following system 

Sv{{k + 1)T) = 6v{kT) +Ta{qv{kT) - eXSv{kT)) 

Y{{k + 1)T) = SiikT) - TaeXhikT) 

maoUsikT) = iB^{q{kT),kT)^)Te^kiqv{kT) 

+ ^k 2 aX{qv{kT) — eXdv{kT))] 

(29) 

where 6 v = [<5i S 2 and S 4 = S 4 — Similarly to the 
state feedback case, if system ( [29] ) is exponentially stable, then 
{q,uj,S) = {q,0,-^q) is an exponentially stable equilibrium 
for (01 (fUji. 











For system (29i, perform the change of variables zi{k) = 
qy{kT), Z 2 {k) = uj{kT)/e, zz(k) = qv{kT) — eX6y{kT), 
Z 4 {k) = S 4 {kT) obtaining 

zi{k) + e’^Z 2 {k) - e^J~^Gi{k,T) 
(B*(fcT)x)T (fcizi(fc) + Ik2a\z4{k)) 


zi{k 

Z2{k 


1 ) 

1 ) 


Z2{k) - eJ-^G2{k,T){B\kTYf 

{k4Zi(k) + Xk2a\z‘i{k)) 
zz^k + l) = zzik) + tT [\z 2 {k) - a\zz{k)) 
-e^J-^Gi{k,T){B^{kT)^y 
[kizi{k) + \k2aXz‘i{k)) 

Z4{k + 1) = Z4{k) — eTaXz4{k) 


The average system of 


defined in ca is given by 


(30) 


zi{k + l) = zi{k) + Z 2 {k) 

Z2{k + l) = Z2{k) - eTJ~'^Lav{T)(kiZi{k)+ k2Z2{k)) 

z:i{k + l) = Z3{k) + eT {kizi{k) + X}^2aXz^{k)) 

Z4{k + 1) = Z4{k) — eT aXz4{k) 

(31) 

Then, a proposition parallel to Proposition holds true, and 
for the average system ( [3T] l the following stability result holds. 

Theorem 2: Assume that fci > 0, /c 2 > 0, a > 0, A > 0 
and Assumption [T] is satisfied. Then there exists T* > 0, and 
for every 0 < T < T*, there exists eg > 0, such that fixed 
0 < T < T* system (311 is exponentially stable for all 0 < 
e < Co- 


Proof: Rewrite in the following compact form 

z{k + 1) = z{k) + eT Ao{T)z(k) (32) 

with 

MT) = 

0 0 o' 

-kiJ-^Lay{T) 0 -\k2aXJ-^Lay{T) 0 

0 \I -aXI 0 

0 0 0 -aA _ 

It will be shown now that = limT->.o X^o{T) is a Hur- 
witz matrix. For that purpose it is useful to introduce the 
continuous-time system w = A^w which in expanded form 
reads as follows 

1 

Wi = -W2 

JW 2 = -L^av (kiWi + \k 2 aXwf) (33) 

Wg = \w2 — aXwg 

W 4 = —aXw 4 


Introduce the Lyapunov function for (33i 


Vsiw) = k4Wj_Ll^Wi 


1 


wIJw 2 + -k2aXwjL°^W3 + -wl 


then 

¥ 3 ( 10 ) = —k 2 a^X^W 3 L°yW 3 — aXwl 


By La Salle’s invariance principle ifThl Theorem 4.4], it follows 
that (33 I is exponentially stable, and thus is Hurwitz. Then, 


the proof can be completed as in the state feedback case. In 
particular, fixing T small enough so that Ao{T) is Hurwitz, 
and letting Pq > 0 such that 


PoAo{t) + AoiryPo = -I 


we have that system ( |3^ is exponentially stable for all 0 < 
e < eg with 


1 

2T\\Ao{TyPoAo{T)\\ 


(34) 


Thus, the following main proposition holds. 

Proposition 3: (Output feedback). Consider the magnet¬ 
ically actuated spacecraft described by ([^ and apply the 
piecewise-constant dynamic output-feedback ( [T^ with k^ > 0, 
A :2 > 0, a > 0, and A > 0. Then, under Assumption [T] 
there exists T* > 0, and for every 0 < T < T*, there exists 
G eo> with eg given by ( [34| , such that fixed 0 < T < T* 
for all 0 < e < ei equilibrium {q,uj, 6) = (q, 0, is locally 
exponentially stable for ([^ ([T4)i. 


V. Case study 

We consider the same case study presented in Q in 
which the spacecraft’s inertia matrix is given by J = 
diag[27 17 25] kg m^; the inclination of the orbit is given 
by incl = 87°, and the orbit’s altitude is 450 km; the 
corresponding orbital period is about 5600 s. Without loss of 
generality the right ascension of the ascending node fl is set 
equal to 0, whereas the initial phase fg (see @) has been 
randomly selected and set equal to fg = 0.94 rad. 

Consider an initial state characterized by attitude equal to 
to the target attitude g(0) = q, and by the following high initial 
angular rate 

w(0) = [0.02 0.02 - 0.03]^ rad/s (35) 


due for example to an impact with an object. 

In Q the continuous-time feedback ( ^ has been designed 
setting fci = 2 10^^, ^2 = 3 10^^, and e = 10“^. Here, we 
keep ki = 2 10^^, k 2 = 8 10^^ and parameters T and e are 
chosen as follows. First, studying numeri cally the behavior of 
the eigenvalues of As{T) with T (see (26i), we determine 
the value T* = 1490 s for which it occurs that As{T) is 
Hurwitz for all 0 < T < T*. Then, the sampling time has 
been set equal to T = 20 to which there corresponds eg = 
1.3 10“^ (see (281). By Proposition]^ it occurs that setting 
e < 1.3 10“^ and small enough, equilibrium {q,uj) = {q, 0) is 
locally exponentially stable for ([^ ([13 1 . Here, proceeding by 
trial and error, we fix e = 10“^. Simulations results reported 
in Fig. [T] show that actually acquisition of the desired attitude 
is achieved. 


A similar approach has been followed to select parameters 
and sampling interval for output-feedback ( [l4| ); the obtained 
results validate Proposition but they are omitted for lack of 
space. 














way. 


0 



orbit 


Fig. 1. Evolutions of {13} with fci = 2 k 2 = 3 10^^, T = 20 s, 

and £ = 10“®. 
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VI. Conclusion 

In this work two magnetic control laws for spacecraft 
attitude stabilizations have been presented. Both magnetic 
control laws possess the feature of generating piecewise- 
constant magnetic dipole moments; the latter aspect makes 
them easier to implement than continuos-time control laws; 
in fact, implementation of continuos-time laws requires a 
discretization process, in which the discretization interval is 
usually selected by trial and error; the latter trial and error 
process in not necessary using the approach here proposed 
since the sampling interval can be determined in a systematic 








































































